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ABSTRACT 

The peaks of 30 optical afterglows and 14 X-ray light-curves display a good anticorrelation of the 
peak flux with the peak epoch: F p oc t p 2 in the optical, F p oc t p 16 in the X-ray, the distributions 
of the peak epochs being consistent with each other. We investigate the ability of two forward-shock 
models for afterglow light-curve peaks - an observer location outside the initial jet aperture and the 
onset of the forward-shock deceleration — to account for those peak correlations. For both models, 
the slope of the F p — t p relation depends only on the slope of the afterglow spectrum. We find that 
only a conical jet seen off-aperture and interacting with a wind-like medium can account for both 
the X-ray peak relation, given the average X-ray spectral slope /3 X = 1.0, and for the larger slope of 
the optical peak relation. However, any conclusion about the origin of the peak flux — peak epoch 
correlation is, at best, tentative, because the current sample of X-ray peaks is too small to allow a 
reliable measurement of the F p — t p relation slope and because more than one mechanism and/or one 
afterglow parameter may be driving that correlation. 

Key words: radiation mechanisms: non-thermal, relativistic processes, shock waves, 
gamma-ray bursts, ISM: jets and outflows 



1 INTRODUCTION 

The power-law fall-off of the GRB afterglow flux is the most- 
often seen feature and was predicted by Meszaros & Rees 
(1997). Long-monitored afterglow light-curves also display 
at least one break. Early (1-10 ks from trigger) breaks in 
the X-ray light-curve have been interpreted as due to energy 
being added to the afterglow blast-wave (Nousek et al 2006) . 
Later (0.3-3 d) breaks in the optical light-curve have been 
predicted by Rhoads (1999) in the framework of a tightly 
collimated afterglow outflow. 

A less-often encountered feature of afterglows is the 
peak displayed by some optical light-curves at early times 
(up to 1 h after trigger). A case by case modelling of those 
peaks could provide a test of the possible models for a light- 
curve peak, as was done, for instance, for the jet model using 
the optical light-curve breaks seen at later times. We have 
investigated models for light-curve peaks only in a general 
sense, by assessing whether a given model can explain the 
strong anticorrelation observed between the peak flux F p 
and the peak epoch t p . Using numerical calculations of the 
reverse and forward-shock synchrotron light-curves for jets 
seen at various angles, we (Panaitescu & Vestrand 2008) 
have found that that model can account qualitatively for 
the F p oc t~ 2 ' 7 peak anticorrelation measured for a dozen 
afterglows with optical light-curve peaks. Analytical results 
for the peak flux and epoch expected at the onset of the 
forward-shock deceleration have been used by us (Panaitescu 
& Vestrand 2011) to test quantitatively the F p — t p relation 
for optical peaks. 



In this work, we continue to investigate the peak mod- 
els mentioned above, by applying analytical results for the 
F p — t p relation expected for the synchrotron emission from 
the forward-shock to afterglows with X-ray peaks. The frac- 
tion of Swift X-ray afterglows displaying a peak is much 
smaller than for optical afterglows because the X-ray peak 
is most often missed, being overshined by the prompt GRB 
emission. We have found about a dozen of X-ray light-curves 
with peaks among the several hundred X-ray afterglow light- 
curves in the XRT database, most of which are accompa- 
nied by a spectral hardness evolution that indicates that the 
emerging X-ray light-curve peak has an origin (= afterglow) 
different than the GRB prompt emission. The importance of 
X-ray peaks lies in that, unlike for optical peaks, the slope 
of the X-ray afterglow spectrum was measured and can be 
used for model testing, as all light-curve peak models lead 
to an Fp oc tp 1 relation with an exponent that depends on 
the spectral slope. As we shall see, the F p — t p peak relation 
is steeper for optical peaks than for X-ray peaks, which we 
will use to further test the two light-curve peak models. 



2 OPTICAL PEAKS 

The sample of afterglows with optical light-curve peaks used 
in this work contains 19 afterglows presented in Panaitescu 
& Vestrand (2011), which are the afterglows with peaks ob- 
served until 2009 and for which the peak fwhm (full width 
at half maximum) is less than a decade (1 dex) in time. We 
add to that set 12 other peaks, some with fwhm larger than 



2 Panaitescu, Vestrand & Wozniak 



1 dex (they were previously considered plateaus or of uncer- 
tain type). For a few afterglows, the peak does not start at 
the first measurement but appears after a flat or a decreas- 
ing optical flux; 110205 displays two clear peaks, both have 
included in the current sample. 

For that set of 31 peaks, the z — 2 optical fluxes 
(F2ev) were calculated assuming an optical spectral slope 
Po = 0.75. Accounting for extinction by dust in the host 
galaxy may lead to intrinsic optical fluxes that are a factor 
up to few brighter, however the 2 = 2 peak fluxes of our 
sample span 5 decades (without 0606014), thus the error in 
the resulting relation between peak flux and peak epoch, 
owing to the (unaccounted for) host extinction, should be 
small. 

To measure the epoch t p and flux F p at the light-curve 
peak, the optical light-curves have been fit with a smoothed 
broken power-law 



F* = f 



1 "V£ 



(1) 



with £ determining the sharpness of the transition between 
the asymptotic F v oc t ar rise and the F v oc t~ ad decay 
(the larger £, the peakier is the light-curve), t„ setting the 
epoch when the transition between these decays occurs (but 
tb is only a very rough estimate of the peak epoch t p ). We 
allow for a shift t for the time when the optical afterglow 
emission begins relative to the GRB trigger (from when t 
is measured), with t > allowing for an optically emitting 
outflow released after the GRB-producing outflow, 

The data for the optical afterglows are shown in Fig- 
ure 1, together with the broken power-law fits. The relevant 
best-fit parameters (a r , a d , t ) to the optical peaks are given 
in Table 1. The peak flux and epoch are not parameters 
of the fit, thus they were calculated from the best- fit pa- 
rameters. With the exception of t , uncertainties were not 
determined, but we provide here an estimate of those uncer- 
tainties for the parameters of interest. The la error of F v 
and t p is less than e = 0.15. The exact value of e(F p ) has 
almost no effect on the following F p — t p relation, while the 
uncertainty on t p is much less than that of t , hence e(t ) 
determines the uncertainty in the true peak epoch (t p — t ). 

The decay index has an absolute error a(ad) < 0.05 
but the rise index error a(a r ) is much larger because to 
its degeneracy with t , both quantities being constrained 
by the measurements during the light-curve rise. The rise 
index a r decreases with increasing t , the uncertainties of 
t and Q r being determined by how well the measurements 
during the flux rise can be described by a concave (holding 
water) "power-law" (for t a < 0) or by a convex (not holding 
water) power- law (for t > 0). For most optical light-curves, 
there are only several measurements during the rise, thus a 
curvature in the power-law rise is often allowed, hence the 
uncertainties of t and a r are often substantial. 

For about 4/5 of the 31 optical peaks, there is no cur- 
vature in the optical flux power-law rise, in the sense that 
the best-fit obtained with t a — (afterglow begins at GRB 
trigger) is statistically as good as the best-fit obtained with 
a free t . More precisely, the relevance of t as a fit pa- 
rameter is established from the F-test probability that the 
increase in the best-fit \ 2 resulting from fixing t — is 



accidental: we consider that t a / is required if the F-test 
probability is less than 10 percent, i.e. there is a higher than 
90 percent chance that the improvement in the best-fit \ 2 
obtained with a free t is real. 

The 1/5 of optical peaks for which the curvature of the 
afterglow rise requires t are listed in Table 2, together 
with the x 2 of the best fits obtained for a free t a , for t = 0, 
and the F-test probability for the t = fit to be acci- 
dentally worse than for the t 7^ fit. There is only one 
afterglow (071010B) for which the optical flux rise is convex 
and requires t < (afterglow begins before the GRB trig- 
ger); for the remaining five cases, t > (afterglow begins 
after GRB trigger) is required by a concave light-curve rise. 

A second argument for the possibility that some opti- 
cal afterglows begun after the GRB trigger is based on the 
model interpretation of afterglow peaks. Numerical calcu- 
lations of afterglow light-curves in the two forward-shock 
models ("off jet-aperture observer" and "deceleration on- 
set") that will be discussed below yield peaks with a fwhm 
of at least 0.6 dex (factor 4 increase in time), depending 
mostly on the post-peak decay index. The effect of a shift 
t > of the afterglow beginning relative to the GRB trigger 
is to increase the fwhm of the peak if time were measured 
from to instead of the GRB onset; thus, all peaks in Table 
1 with fwhm < 0.6 would require t > if they are to be 
accounted for the peak models considered below. There are 
five such narrow optical peaks in Table 1 and, for two of 
them (see Table 2), we do find that t > is also required 
by the best-fit to their rises. 

Figure 2 shows the broken power-law fits to the 31 op- 
tical peaks and the best-fit power-law relation between peak 
time t p and peak flux F p , obtained by calculating the coef- 
ficients a and b that minimize 



^2 _ W°sFp,i - (a + blogt p ^)] 



4^ <7 2 (logF p ,;) +& 2 <7 2 (log t p ,i) 



(2) 



For t = (i.e. the peaks labeled in Fig 2 by the GRB 
name), the best-fit is F p oc tp 2 ' 15±0 ' 07 , assuming the same 
relative error e(t p ) = e(F p ) = 0.15 for all afterglows. The 
exponent changes by about 0.10 if either e(t p ) or e(F p ) is 
varied by a plausible 0.05, thus the full la uncertainty of 
the t p exponent is about 0.12 (adding in quadrature the two 
uncertainties). 

This best-fit is significantly less steep than the one 
previously reported by us (Panaitescu & Vestrand 2011), 
F p oc t p 3 ' 2±0 ' 2 , for a set of 16 optical peaks, owing to the in- 
clusion of a few afterglows which peak at later epochs. It is, 
however, consistent with the correlation found by Liang et 
al (2012), Fp oc t~ 19±0 - 3 , for a set of 39 afterglows. Twenty- 
three of the peaks in Table 1 are in Liang's sample, for the 
other 16 afterglows either the redshift is not known or the 
optical measurements do not allow a clear identification of 
a peak in the light-curve. 

The linear (Pearson) correlation coefficient 



r(x,y) = 



< xy > — < x >< y > 



(3) 



for the optical peaks shown in Figure 2 is r(log F p , log t p ) = 
—0.77 ± 0.01. For 31 points, the corresponding probability 
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Figure 1. Optical measurements for 30 afterglows and smoothly broken power-law fits to 31 peaks, obtained with the function given in 
equation (1). Note that not all existing measurements have been fit, but only those that represent a light-curve peak. 
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Table 1. Parameters of the smoothed broken power-law best-fit (eq 1) to z = 2 optical light-curves with peaks (Fig 1) and other 
quantities that characterize those peaks 
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(1) number of optical measurements, (2) epoch of first measurement, (3) power-law rise index of the best-fit for t > (uncertain 
because of its degeneracy with the time-origin t ), (4) epoch of light-curve peak (is usually close to the break-time t b ), (5) width 
of the peak at half maximum, in log scale, and for t a = (fwhm is a substitute for the smoothness parameter £) (6) 2eV flux at 
light-curve peak (is usually within a factor 2 to the normalization factor /), (7) best- fit power-law decay index, (8) epoch of last 
measurement used for fit, (9) reduced \ 2 of the best-fit obtained with t > 0, (10) references for data. 

Table 2. Optical afterglows for which t ^ provides a better fit to the light-curve rise than t a = 0, i.e. for which the afterglow beginning 
is not the GRB trigger 
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(1) number of measurements during the afterglow rise, (2) reduced 
X 2 of the best-fit obtained with free t , (3) lower limit of the 90 
percent confidence level (cl) on t a , (4) upper limit of the 90 percent 
cl on t , (5) reduced \ 2 °f tne best-fit obtained with t = 0, (6) 
chance probability (in percent) that t a = 0. 
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Figure 2. Smoothly broken power-law fits to 31 optical peaks. 
The peak location with time measured since the GRB trigger is 
labelled by the GRB name; the best-fit in log-log scale is shown 
with a red line. Triangles show the first measurement for sev- 
eral optical afterglows whose rises were missed, i.e. those peaks 
occurred before the first measurement. 

that this correlation occurred by chance (i.e. in the null hy- 
pothesis, that the peak epoch is not correlated with the peak 
flux) is logp = —6.5. 

If we allow for a time-shift r = \\to + ti" 111 ^] be- 
tween the GRB trigger and the onset of the optical after- 
glow, with [ti"" n ^ — t max ^] the 90 percent confidence level 
(cl) on the afterglow onset epoch, determined by a variation 
of A^ 2 = 2.7 around the best-fit obtained with t = 0, then 
r[log F p , log(£ p — t)] = —0.80 ± 0.02 and the corresponding 
chance probability logp = —6.9 is 3 times lower than for the 
F p — t p correlation. The new best-fit is 

F { p° pt) = 3.6 x 10 s {tp - r) -7 " (mJy) , 7o = 2.00 ± 0.08 (4) 

assuming a relative error e(F p ) — 0.15 for all afterglows 
and using the uncertainty of peak-time shift, <t(t) = 
0.60(ti max) - t ( min) )/2, as the la error for t p - r. These 
results are almost unchanged if one uses t of the best-fit 
instead of the middle of the 90 percent cl. 



3 X-RAY PEAKS 

X-ray afterglows light-curves displaying a peak are rarely ob- 
served by Swift-XRT, in the sense that we found only such 



14 cases among the several hundred of X-ray light-curves 
in the UK Science Data Center repository. One reason for 
that rarity is that the prompt emission and its tail over- 
shine the afterglow emission until tens-hundreds of seconds 
after trigger. Another reason is that we have retained only 
afterglows whose emergence after the GRB tail was suffi- 
ciently well sampled. A strong spectral hardening is usually 
observed at the emergence of the X-ray afterglow, but such a 
spectral evolution was not a criterion for selecting the X-ray 
afterglows. 

Some of the X-ray light-curve peaks are displayed by 
the 0.3-10 keV band light-curves found at the Swift-XRT 
light-curve repository (Evans et al 2007, 2009, 2010), but 
others display a peak only in the XRT processed 10 keV 
light-curve (for all, we have calculated the 1 keV flux using 
the X-ray spectral slope reported at the Swift-XRT spec- 
trum repository - Evans et al 2007). The reason for which 
some 0.3-10 keV light-curves do not display a peak although 
the processed 10 keV light-curves do so is that, just after the 
emergence of the afterglow rise, the softer GRB tail emis- 
sion contributes more to the 0.3-10 keV flux than to the 10 
keV emission, relative to the harder afterglow flux. Thus, for 
afterglows peaking shortly after their emergence from under 
the GRB tail, the rise of the 0.3-10 keV afterglow flux can be 
masked by the GRB tail. Whenever they displayed a clear 
peak, we used the 0.3-10 keV light-curves instead of the 10 
keV fluxes, because the flux errors of the 10 keV processed 
light-curves are underestimated, as the errors of the inter- 
polated hardness ratios (or spectral slopes) were not taken 
into account in the calculation of the 10 keV flux from the 
0.3-10 keV count rate. 

The 14 X-ray afterglows with peaks are shown in Figure 
3 and some of the parameters obtained with the smoothly 
broken power-law fit are listed in Table 3. The allowed range 
for the X-ray emission onset epoch t were calculated only 
for t > 0, the true 90 percent cl being even larger, owing to 
the (ct r ,t ) fitting degeneracy. For most of X-ray afterglows, 
the 90 percent cl on t extends from trigger to near the 
first measurement, with only two exceptions: for 080319C, 
t > 220 s, while for 091029, t > 410 s, the probability 
(calculated from the corresponding \ 2 ) that t = being 
around 10 percent for both afterglows. 

The 90 percent cZs on t are often quite large, being in 
some cases the entire interval [0, ti] from GRB trigger and 
until first measurement, because of the degeneracy between 
t and the rise slope a r . Narrower ranges for t are obtained 
if a r is restricted to, for instance, the range [2, 3] spanned 
by the expectations for the pre-deceleration forward-shock 
model (and for homogeneous ejecta). However, because some 
peaks display a rise slower than t 2 if time is measured from 
GRB trigger, it is necessary to allow t < 0, so that the rise 
of such afterglows becomes faster when time is measured 
from t . With the time shift r determined as above, but 
using the 90 percent els on t calculated from fits with the 
2 < a r < 3 restriction, the F p — t p correlation is significantly 
stronger (logp = -7.7) and the best-fit is F p oc ^ 1 - 90±0 07 . 

The last column of Table 3 lists the slope of the after- 
glow X-ray spectrum reported at the Swift-XRT spectrum 
repository. The error-weighted average of those 14 slopes is 

(F v oc 1.04 ±0.02 (5) 
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Figure 3. Light-curves of 14 X-ray afterglows with peaks and broken power-law fits to those peaks. 



Table 3. Parameters of the smoothed broken power-law best-fit to 14 X-ray light-curves with peaks (Fig 3), at z = 2, and other quantities 
of interest describing those light-curves, as defined in Table 1. Last column gives the slope of the X-ray spectrum F v oc u~ (and its 90 
percent cl). 
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fwhm 
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At 
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(s) 




(fl) 


(dex) 






(ks) 






050824 


10 


11000 


0.9 


47000 


1.1 


0.019 


4.47 


540 


0.64 


0.95(.18) 


051016B 


46 


390 


0.8 


1500 


0.9 


0.32 


0.74 


60 


0.85 


0.89(.13) 


060926 


21 


185 


4.0 


310 


0.6 


17 


1.36 


36 


0.90 


1.00(.29) 


070103 


27 


200 


3.6 


540 


1.0 


11 


1.45 


106 


0.87 


1.31(.24) 


080310 


12 


740 


1.1 


2000 


0.9 


1.5 


0.75 


30 


3.41 


1.09(.07) 


080319C 


58 


260 


1.8 


350 


0.6 


430 


1.82 


350 


1.08 


0.94(.32) 


081028 


19 


7800 


0.4 


14600 


0.6 


3.3 


2.31 


400 


0.94 


1.03(.07) 


090205 


31 


230 


1.1 


490 


0.8 


9.3 


0.78 


13 


0.80 


1.03(.15) 


090429B 


20 


46 


1.3 


180 


1.0 


75 


1.44 


39 


1.08 


1.00(.24) 


091029 


53 


830 


0.2 


1500 


1.0 


3.4 


1.34 


1300 


0.98 


1.09(.10) 


100418 


26 


1070 


0.7 


76000 


1.6 


0.013 


3.5 


3700 


1.59 


1.27(.40) 


100513 


22 


360 


0.2 


1900 


1.1 


6.1 


0.96 


170 


1.11 


1.27(.27) 


100901 


42 


12000 


0.7 


34300 


0.8 


0.87 


1.49 


1640 


1.63 


1.09(.07) 


110213 


20 


160 


0.9 


1300 


1.0 


44 


2.09 


360 


1.66 


0.99(.07) 



error being la, and with a xt — 1-13 for all slopes being 
consistent with the average fi x given above. 

Figure 4 shows the broken power-law fits to the 14 X- 
ray peaks together with the best-fit to their peak fluxes and 
epochs. For t — o, the linear correlation coefficient of the 
peak fluxes and peak epochs is r(\ogF p ,logt p ) — —0.83 ± 



0.01, chance probability logp = —3.7, and the best-fit to 
the peak locations is F p oc t~ 1,59±0 ' 07 . For the peak epochs 
shifted by the center of the 90 percent cl on t > 0, we obtain 
r[logF p ,log(tp - r)] = -0.85 ± 0.02, logp = -3.9 and 

F ( X ray) = ^ fe _ r) - 7x (mJy) _ ^ = Lg g ± Q Qg (g) 




time (s) for z=2 

Figure 4. Broken power-law fits to the peaks of 14 X-ray after- 
glows and the power-law fit to their peak fluxes and epochs (blue 
line). Peak location is indicated by the afterglow name. 



which is shallower than for optical light-curves (eq 4): 

7* - To = -0.37 ± 0.12 . (7) 

To compare the distributions of optical and X-ray peak 
epochs, one can calculate 



2 _ x - (y/N x /N m, - ^N /N x m, x ) 2 



(8) 



with m the number of optical or X-ray peaks in time-bin 
i and N the total number of optical or X-ray afterglows. 
For our two samples, the resulting \ 2 = 6.7 for 7 dof cor- 
responds to a 50 percent probability that the optical and 
X-ray peaks shown in Figure 5 are drawn from the same 
distribution. 

ORIGIN OF PEAK FLUX - 
PEAK EPOCH ANTICORRELATION 

As indicated in Figure 2, for some optical afterglows, 
we have seen a decaying flux since first measurement, thus 
the peak of the optical light-curve occurred at an earlier 
time than for the afterglows listed in Table 1. This means 
that equation (4) represents only the bright/late edge of the 
distribution of all optical light-curve peaks in the (F p ,t p ) 
plane. 

There are two mechanisms that can produce a peak 




10 10 

time (s) for z=2 

Figure 5. 31 optical peaks (red symbols) and 14 X-ray peaks 
(blue) with the peak epoch measured from the middle value of 
the 90 percent confidence level on f , determined by fitting light- 
curves with the function in equation (1). The histograms at the 
bottom show the distributions of optical (red) and X-ray (blue) 
peak epochs. 



in the afterglow light-curve: 1) an observer location that 
is (for some time) outside the aperture of the relativistic 
afterglow jet, such that the received afterglow flux rises when 
the observer is still outside the 1 /F aperture of the afterglow 
emission, T being the ever-decreasing jet Lorentz factor, and 
2) a dynamics of the afterglow-producing shock such that 
the flux of that shock's emission rises for some time. 

In the former case, the peak of the afterglow light-curve 
occurs when, owing to the gradual increase of the aper- 
ture of the Doppler-boosted emission cone, the center of the 
jet becomes visible to the observer. In the latter case, the af- 
terglow rises as more energy is added to the afterglow shock, 
the light-curve peak corresponding to a change in the rate 
at which that energy is injected into the shock. A sudden 
change in the injected power occurs naturally if the ejecta 
are contained in a well-defined shell, so that the dynamics of 
both shocked fluids (ejecta and swept-up medium) changes 
(i.e. deceleration sets in) when the reverse shock crosses the 
entire ejecta shell. 

For either model for light-curve peaks (off-jet observer 
location or onset of deceleration), we will search for the pa- 
rameter whose variation among afterglows can explain the 
existence of the bright /late edge of peak fluxes and epochs in 
the (F p , t p ) plane (i.e. can induce the F p — t p correlations of 
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eqs 4 and 6), and that can yield a steeper F p — t p dependence 
for optical peaks than for X-ray peaks. 



4 OFF-JET OBSERVER LOCATION 

A structured or dual outflow model was first proposed by 
Berger et al (2003) to explain the radio emission of GRB af- 
terglow 030329. The existence of many afterglows with de- 
coupled optical and X-ray afterglow light-curves (i.e. with 
different decay rates at these two frequencies or with chro- 
matic X-ray breaks) has revived interest in this model (e.g. 
Racusin et al 2008). 

The basic feature of this model for light-curve peaks 
is that the observer is located outside the opening of the 
afterglow jet, 9j, at an angle 9 > 9j. That is a somewhat 
unappealing feature, as it implies that the GRB emission 
is produced by a different jet, whose aperture includes the 
direction toward the observer, allowing us to localize the 
burst and follow its afterglow. A more palatable version of 
it is an outflow with an non-uniform angular structure (i.e. 
distribution of ejecta energy with direction in the jet), such 
as a bright core moving toward the observer and emitting 7- 
rays during the prompt phase and an envelope surrounding 
the core an dominating the afterglow emission at later times. 

For simplicity of derivations, we focus below on the off- 
aperture jet, but the results for the peak flux and epoch 
should be the same for an envelope outflow. The emission re- 
ceived by the observer is that in the comoving frame, boosted 
relativistically, F v = T> 3 F^ T>: where T> is the relativistic 
boost factor 



V - 



2F _ \2/{V6 2 ) (T>8- 1 ) 



r(i-vcos0 o ) ~ r 2 <9 2 + i ~ \2r (r < e- 1 ) 



(9) 



r being the jet Lorentz factor, and assuming F 1 and 
do <C 1. At early times (i.e. when T > 6C 1 ), the Doppler 
boost increases while at later time (i.e. when F < 0~ 1 ), the 
Doppler factor decreases, both owing to the jet deceleration. 
The behaviour of that factor is convolved with the decrease 
of the afterglow flux in the comoving frame, at frequencies 
above the peak of the comoving-frame (synchrotron) emis- 
sion spectrum, also due to the continuous jet deceleration. 
For r > 0"\ the increase of T> dominates and the observer 
sees an increasing afterglow flux, until F = O^ 1 , when the 
afterglow light-curve peaks, being followed by a decreas- 
ing flux when F < Thus the afterglow peak occurs 
when the ever-widening r _1 -opening cone of relativistically- 
boosted emission contains the direction toward the observer, 
at which time the Doppler boost is 



V{t p ) = F{t p ) = 9- 1 . 



(10) 



The observer offset angle 8 is the only parameter to 
which makes sense to attribute the the F v — t p correlation, 
as an increasing angle 9 yields a later peak epoch t p and 
a lower peak flux F p (see below). To calculate the F p — t p 
relation induced by a variable angle O , one has to relate F p 
and t p to the offset angle. The observer-frame arrival time 
t of the photons emitted by a relativistic source moving at 
Lorentz factor F(r) and angle 9 a relative to the direction 
toward the observer is found by integrating 



cdt = — — dr cos t 

v 



^(r- 2 + 9>) 



(ii) 



To calculate the afterglow flux at t p , we restrict our at- 
tention to the forward-shock synchrotron emission (i.e. the 
emission from the shocked ambient medium), and leave the 
derivation of the F p — t p relation for the reverse-shock emis- 
sion to the fans of that model. In the comoving frame, the 
forward-shock synchrotron emission peaks at frequency 



/ 2 D r 3 1/2 

i>i oc 7, B oc 1 n ' 



(12) 



with 7i oc T being the typical post-shock electron energy, 
B 2 oc Fn' oc T 2 n the magnetic field energy density, n' oc Fn 
the post-shock particle density, and n the density of the 
circumburst medium. The comoving frame synchrotron flux 
density at v[ is 



F'i = F'<x N e B oc MFn 1 ' 2 



(13) 



where N e oc M is the number of forward-shock electrons 
and M the mass of the energized ambient medium. The 
synchrotron spectrum has a "cooling" break at the charac- 
teristic synchrotron frequency v' c of the electrons whose ra- 
diative cooling timescale (t' rad = E' el /P!,y n oc jc/^B 2 ) = 
7 < T 1 .B~ 2 ) equals the dynamical timescale (t' dyn oc r/T), thus 
7 C oc FB^r' 1 and 



v c oc ~j 2 B oc F 2 B 3 r 2 oc F 1 n 



1^-3/2^,-2 



(14) 



In the observer frame, the characteristics of the syn- 
chrotron spectrum are 

F = V s F' , Ui = Vvl , v c = Vv' c (15) 
and the afterglow flux at an observing frequency v is 

F v {y c <v<v i ) = F (^) oc V 7 ' 2 F'Mf 12 (16) 
F v (yi <v<v c ) = F oc V^F'Mf (17) 

F v (i>i, v c <v)=f q ^ i/2 ^y x v^F^r^y\i8) 

where /3 is the slope of the afterglow spectrum at frequency 
observing v. The case given in equation (18) applies only 
to optical light-curves because the implied spectrum (F v oc 
v^ 1 / 2 ) could be compatible with that of optical afterglows 
but is clearly harder than observed for X-ray afterglows (eq 
5). 



4.1 Synchrotron forward-shock emission from 
conical jets 

For an adiabatic jet that does not spread laterally (because, 
for instance, the jet is embedded in a confining envelope), 
conservation of energy during the interaction between the 
jet and the ambient medium reads E oc F 2 M with E the jet 
energy and M oc n(r)r 3 oc r 3 ~ s . Then, for a n oc r~ s (with 
s < 3) radial distribution of the ambient medium density, 
the jet dynamics is given by 



F oc r 



-(3-s)/2 



(19) 
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and equation (11) yields 

t = I( i l 7 r- a + fl2) (20) 

hence the epoch of the light-curve peak, t p = t\ r _„-i, sat- 

isfies t p oc r p 6l, where r p = r\ T=0 -i oc is the jet 

radius at which the light-curve peaks, the last relation fol- 
lowing from equation (19). Thus 



t v oc I 



-2*)/(3-<>) 



(21) 



Eliminating o between the expressions for r p and t p , it fol- 
lows that 



r„ oc t 



(22) 



Substituting n(r p ), M(r p ), and T(r p ) in equations (12) 
- (14) and using equation (22), we get 

K(t P ) « r p /2_S > "i(* P ) « r^ 9/2 , i»c(tp) oc r p ~ 1/2 . (23) 

The dependence of the afterglow peak flux F p = F v (t p ) on 
the light-curve peak epoch t p can be now calculated from 



(s-3)/(8-2s) 



and us- 



equations (15)-(18), with T> oc 0„ cc tp 
ing the spectral characteristics given in equation (23), with 
r p (t p ) from equation (22). 



4-1.1 Homogeneous medium (s=0) 

With the above-mentioned substitutions, one obtains the 
following exponents for the F p oc t p 7 correlation for after- 
glow light-curves from off-aperture jets interacting with a 
homogeneous medium (s = 0): fe(v c < v < Vi) = 1, which 
is too small compared with the values measured for optical 
and X-ray peaks, and 

3 1 

Vi < v < v c : 79 = -(2/3+1) , v i: v c <v: f 9 = ^(6/3+1) .(24) 

For the former case, the average measured X-ray spectral 
slope (eq 5) implies that f x = 2.3, which is too large than 
what is measured for X-ray peaks (eq 6). The second case 
leads to f x = 1.81 ± 0.03, which is almost compatible with 
observations. However, this case cannot account for the slope 
7o measured for X-ray peaks, if the optical and X-ray af- 
terglow emissions arose from the same outflow (same elec- 
tron population): for Vi,v c < v Q (i.e. j3 = /3 X ), we ex- 
pect 7o = fx, which is inconsistent with observations; for 
Vi < v < v c < v x , (i.e. p x — j3 + 1/2), equation (24) im- 
plies that f x — 7 = \.h[fix ~ /So) — 0.5 = 1/4, contrary to 
observations (eq 7). 



4-1.2 Wind-like medium (s=2) 

Similar to above, for a medium having the n oc r~ 2 radial 
stratification expected for the winds of massive stars (as the 
progenitors of long GRBs), we obtain fe(v c < v < Vi) — 
0.75, which is smaller than measured for optical and X-ray 
peaks, and 

3 3 
Vi < v < v c : fg = -f) + 1 , Vi, v c < v : fe = -ft . (25) 



In the former case, the average X-ray spectral slope leads 
to fx — 2.6, which is larger than observed, while the latter 
case leads to f x = 1.56±0.03, which is consistent with obser- 
vations of X-ray peaks, but too small for the optical peaks 
slope. However, if Vi < v < v c < v x , then equation (25) 
yields 7^—70 = 1.5(fi x —j3 ) — 1 = —1/4, which is marginally 
consistent with the measured value (eq 7) . Thus, a jet seen 
off aperture, interacting with a wind-like medium, can ac- 
count for both the optical and X-ray peak slopes provided 
that the synchrotron cooling frequency is between optical 
and X-rays. 



4-1.3 Jet half-opening and observer offset angle 

In deriving the above scalings for the flux, we have ignored 
any intrinsic afterglow parameter and kept only the depen- 
dencies on the observer's offset angle O because we are in- 
terested in the F p — t p relation induced by the variation of O . 
The ignored afterglow parameters, kinetic energy per solid 
angle E and ambient density n, are not the same for all af- 
terglows, and their variation of E and n among afterglows 
yields a scatter in the (F p ,t p ) the plane, along the relation 
induced by the variation of o . 

Equation (21) implies that, aside from that scatter in- 
duced by the variation of E and n, the offset angle 9 must 
span a range of 1 dex, to account for the 3 dex spread in 
observed peak time t p (Fig 5). For typical afterglow ener- 
gies and ambient medium densities, the jet Lorentz factor 
at the earlier peak epochs shown in Figure 5, of about 100 
s, is expected to be T ~ 100, which implies that the small- 
est offset is 0.5 deg. To obtain a light-curve peak, 
the jet opening must be smaller than the observer's offset, 
thus the earliest peaks correspond to the narrowest jets, with 
Q(mm) ^ q g j£ a jj j e ^ g game p en i n g ) then the 

largest offset angle would be 0^ o max ^ 5 deg. However, the 
latest occurring peaks are not necessarily from narrow jets 
seen at a larger offset angle, they could also be from wider 
jets: at 10 ks, the jet Lorentz factor is F ~ 10, hence the 
offset angle is o ~ 5 deg and the widest jet could have 
g(m,a X ) ^ g Therefore, the range of peak epochs shown 
in Figure 5 suggests that afterglow jets have half-openings 9j 
from 1/2 to 5 degrees and are seen at an angle o = (1+N)6j, 
with N < 10. 



4.2 Synchrotron forward-shock emission from a 

spreading jet interacting with a homogeneous 
medium 

If the afterglow jet spreads laterally unimpeded, then the 
jet dynamics changes when the lateral spreading (at the 
sound speed) increases significantly the jet opening, lead- 
ing to other dependencies of the peak flux F p on the peak 
epoch t p than derived above for a conical jet. Defining rj to 
be the radius at which the jet Lorentz factor has decreased 
to Fj = 6~ , where 0j is the initial jet half-angle, the jet 
dynamics is given by (Panaitescu & Vestrand 2011) 



T(r) = 



-(— -1) 



3/2 (r < r 3 ) 



(26) 
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for a homogeneous medium (the jet dynamical equations 
cannot be solved analytically for a wind- like medium). At 
r < Tj, the lateral spreading is negligible and the jet dynam- 
ics is the same as for a conical jet but at r > rj the lateral 
expansion of the jet more than doubled its initial aperture 
and an exponential deceleration sets in (Rhoads 1999). 

With the exception of observer locations just outside 
the jet opening, the afterglow light-curve peak occurs after 
rj (i.e. r v > rj), when the jet lateral spreading is significant. 
Equation (11) can be integrated first over the power-law 
jet deceleration and then over the exponential deceleration, 
yielding 



t = 



4r 2 



(27) 



For do S> 0j, the jet becomes visible to the observer at radius 
r p > rj, thus 



(28) 



where we used the definition of the afterglow peak time: 
r(t p ) = 6~ 1 . That condition and the dynamics of the jet at 
r > rj (eq 26), imply that 



(29) 



This result can be simplified if it is assumed that all jets 
have the same 9j and rj, which is equivalent to assuming 
that all jets have the same 8j and ratio E/n. This simplifi- 
cation is motivated by that we are searching for the F p — t v 
dependence arising only from varying the observer location, 
and not from varying the jet dynamical parameters (which 
is the subject of the next section, but in an other model for 
light- curve peaks) . Thus, 



r p e r *> /r i oc O oc (t p /r p ) 1/2 



(30) 



the last scaling resulting from equation (28). 

It can be shown that, for a homogeneous medium, M oc 
r 2 exp(2r/rj) at r > rj. Substituting it and the jet dynamics 
(eq 26) in equations (12) - (14), we get 



F' {t p ) oc r p e r * /r i , 
v[{t p ) oc r- 3 e- 3r ^ 
v' c (t p ) oc r p 1 e rp/r * . 



(31) 
(32) 
(33) 



From here, one can calculate F p = F v (t p ) by using equations 
(16) - (18), with T) oc O^ 1 and by substituting exp(r p /rj) 
with the aid of equation (30): 



F p (y c < v < m) oc t p 1 



FJ Vi <v<v c )ct rf+X^+V 



F p (yi, v c < v) oc r 2 p 1 t p 2/3 



(34) 
(35) 
(36) 



Ignoring the r p factors (because r p has a weak, sub- 
logarithmic dependence on t p , according to eq 30), one ob- 
tains: 7e(^ c < v < Pi) = 1, which is inconsistent with the 



F p — t p correlation measured for optical and X-ray peaks, 
and 



Vi < v < v c : je — 2/3 + 1 , Vi , v c < v : je — 2/3 



(37) 



In the former case, the average X-ray spectral slope leads to 
an exponent 7 = 3.1, inconsistent with that measured for 
X-ray peaks, while the latter case yields 7 = 2.08 ± 0.04, 
which is too large for X-ray peaks but consistent with that 
of optical peaks, thus this case can explain the the F v — t v 
relation for optical peaks if optical were above the cooling 
frequency. If the cooling break were always between optical 
and X-rays, then the spreading-jet model expectation is -y x — 
7o = 2(/3 x — fio) — 1 = 0, inconsistent with observations. 



5 ONSET OF DECELERATION 

In this model, for afterglow peak occurs when the reverse 
shock has finished crossing the ejecta shell and the mo- 
tion of the shocked fluids starts to decelerate (or to decel- 
erate faster than before). We focus only on the forward- 
shock emission, which depends on fewer afterglow parame- 
ters than the reverse-shock's. Depending on the initial geo- 
metrical thickness A of the ejecta shell, the observer-frame 
deceleration timescale (i.e. the afterglow peak epoch) is ei- 
ther A/c (thick ejecta shell/relativistic reverse-shock case) 
or is determined by the parameter set (E, Fo, n) (thin ejecta 
shell/semi-relativistic reverse-shock), where E is the ejecta 
initial kinetic energy per solid angle and To the ejecta initial 
Lorentz factor. 



5.1 Forward-shock synchrotron emission for a 
semi-relativistic reverse-shock 

If the ejecta is shell is sufficiently thin (or, equivalently, suf- 
ficiently dense), the reverse-shock is semi-relativistic, the 
shocked fluids (ejecta and ambient medium) move at a 
Lorentz factor T lower than but close to that of the un- 
shocked ejecta To, and the time it takes the reverse shock 
to cross the ejecta shell is close to the time it takes the 
forward-shock to sweep-up a mass of ambient medium equal 
to a fraction T^ 1 of the ejecta mass. 

Thus, the deceleration radius is defined by M(rd) = 
E/r' 2 ,, from where 



ra oc 



E 

fcrf 



l/(3-s) 



for an external medium of proton density n(r) = kr 
the observer-frame deceleration timescale is 



ta = 



2 C rg 



(38) 



and 



(39) 



The condition for a semi-relativistic shock is simply A < ct d , 
which can be better expressed as an upper limit on To. 

The flux at ta can be calculated starting from equations 
(12) - (15), with two modifications. First, the Doppler factor 
is just the Lorentz factor of the shocked fluid at td, which is 
~ r () . Second, F — YqF' because, owing to the relativistic 
beaming, the observer receives emission from a region of 
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opening 9 = T 1 that contains a fraction T 2 of the number 
of electrons iV e oc r 3 n for a spherical source. Therefore 



E 



E 



r p OC fcr g , tp OC ^4 



F 



r -F„ cc iV e Br oc Mron 1/2 



oc r n 



2 3/2 3 



-p / -1^4 1/2 / -3/2 -2 

= f o^i oc 1 n ' , v c = 1 o^ c oc n ' r 



5.i.-Z Homogeneous medium 

For n(r) = no, equations (38) and (39) lead to 

1/3 / ^ \ 1 / 3 



(40) 
(41) 



-<- :t " »i(t p ) oc B-^fc 3 / 2 , 

l-/3p6/3-2 fc (3/3-l)/2 
1/2 F „ 1.1/4 



(49) 

2 1.-5/2, 



Fo(t p ) oc rgr^, i/i(t p ) oc s^rgr^, i> c (* p ) oc £tvat°^(50) 

(51) 
(52) 



f p (^,^ c < !/) oc s'-T^fc 



i B oc 



E 



(42) 



By substituting r = r d with r d from equation (38) in equa- 
tions (40) and (41), the spectral characteristics at the peak 
epoch t p = td satisfy 

F (t p ) oc En /2 , Vi(t p ) oc T^nJ /2 , v c (t p ) oc £- 2/3 r 4/3 n ( 7 5/6 .(43) 

Then, for the three orderings of Vi, v, and v c that yield a 
peak for the afterglow light-curve at the onset of deceleration 
(i.e. a rising light-curve before t d and a falling-off flux after 
that), we get the peak flux 

F v {vi < v < v c ) oc F vf\ tp oc ETf n (p+1)/2 (44) 
Fp{vi, v c < v) oc J>yV-°- 5 | tp oc £ 2 / 3 r^- 1 / 3 »n^- 1 '/ 6 (45) 
F p (i/ C < v < Vi) oc J> c 1/2 | tp oc S 2 / 3 r 2/3 n 1/12 . (46) 

Comparing with equation (42), the above results show 
that a variation of E among afterglows induces a F p — t p pos- 
itive correlation (unlike the observed anticorrelation) . For 
the case in equation (44), varying n and To produces anti- 
correlations F p oc 7 with 



3 3 
Vi < v < v c : 7„ = -(/? + 1) , 7r = -/3 



(47) 



thus the average X-ray spectral slope leads to indices 7 x>n = 
3.1 and 7^ = 1-56 ± 0.03, respectively. The latter case is 
consistent with the slope measured for X-ray peaks, however, 
in that case, /3 = p x and 7 D = j x , thus it cannot account 
for the measured optical peaks slope. The measured spread 
of 3 dex in peak epochs would require that To ranges over 
about 1 dex (from eq 42) . 

For the case given in equation (45), varying n and To 
leads to the same anticorrelation slope 



v c < v : 7„ = 7r = -(3/3 - 1) 



(48) 



thus the measured X-ray spectral slope implies -y x = 1.1, 
which is too small. 

For the case in equation (46), the corresponding after- 
glow spectral slope j3 — 1/2 is harder than usually measured 
in the X-ray and the resulting indices 7„ = 7r = 1/4 arc 
well below those measured. 



5.1.2 Wind-like medium 

For n(r) = kr~ 2 , the equations above become 



F p (v c < v < Vi) oc E'^Fok 1 



From equation (51), the expected exponents of the F p — 
tp relation are 



Vi,v c <v : 7is = - 1 , 7 fc = 7r = - (3/3 - 1) (53) 

thus the measured average X-ray spectral slope implies 
7x,b = and 7 I . fc / r = 1-1, respectively, both being less 
than measured. For the case given in equation (52), 7^ = 
7r = 1/4, neither of which is compatible with observations. 
For Vi < v < v c the onset of deceleration does not yield 
a light-curve peak, the afterglow flux decreasing even the 
reverse-shock energizes the ejecta.* 



5.2 Forward-shock synchrotron emission for a 
relativistic reverse shock 

If the ejecta shell is sufficiently thick, the reverse-shock is 
relativistic and the shocked fluids move at a Lorentz factor 
r that can be well below the Fo of the incoming ejecta. It 
can be shown (Panaitescu & Kumar 2004) that the reverse- 
shock reaches the trailing edge of the ejecta shell at a radius 



Td oc 



/£A\ 1/(4 - s) 



V k J 

when the radiating fluid moves at 

E ,1/(8-23) 



T oc 



(fcA 



(54) 



(55) 



corresponding to an observer-frame time td = rtj/(2cr 2 ) ~ 
A/c, i.e. the (new) observer-frame deceleration timescale is 
about the lab- frame geometrical thickness of the ejecta shell. 
The spectral characteristics at ra are those of equation (40) 
and (41) with Y of equation (55) instead of Fo. Given that 
t d depends only on A, to derive the peak flux - peak epoch 
relation requires only the dependence of the spectral char- 
acteristics on A. 



For a homogeneous medium, 



r d oc A 1/4 , r oc A- 3/s 



F {tp) oc A , Vi(t p ) oc A 3/2 , Vcitp) oc A 1/2 



(56) 
(57) 



and the slopes of the F p — t p relation induced by the variation 
of A among afterglows are 



Vi < v < v c 



7A = -/ 



< v < Vi : 7a = 1/4 



(58) 



* A light-curve peak would be produced only by the synchrotron 
peak frequency Vi falling through the observing band; however 
that is a different light-curve peak mechanism and implies a hard 
F v oc i/ 1 / 3 X-ray spectrum before the peak, which is not seen by 
XRT for any of the afterglows in Table 3. 
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For the first case, the expected slope is y x = 1.56 ± 0.03, 
which is consistent with observations, but it implies that 
7o = 7x, hence it does not provide an explanation for the 
larger exponent measured for optical peaks. For Ui, v c < v, 
the onset of the deceleration does not yield a light-curve 
peak, as the afterglow flux decreases while the reverse-shock 
crosses the ejecta. 

The same issue (deceleration onset does not yield a 
light-curve peak) exists for a wind-like medium, whatever 
is the fi x > spectral regime (as measured by XRT). 



6 CONCLUSIONS 

We have identified 31 optical peaks (Fig 2) with measured 
redshifts and we fit them with a smoothly broken power-law, 
to determine their peak times and epochs more accurately, 
and also to assess if a time-shift t between the beginning of 
the afterglow emission and the GRB trigger is required. We 
found good evidence that t < for only one afterglow and 
t > for five other, while for the rest t =fi does not pro- 
vide a statistically significant better fit. Because we fit the 
afterglow light-curves with power-laws, finding that the af- 
terglow beginning is different than the GRB trigger happens 
mostly when the afterglow rise (with time measured from 
burst beginning) is curved in log-log space. Consequently, 
such determinations of t rely on the untested assumption 
that the intrinsic afterglow rise is a power-law. ' 

The fc-corrected optical light-curves manifest a strong 
correlation between the peak flux F v and the peak epoch 
t p . The best fit to the afterglow peaks in the log- log plane is 
F p cc tp 216 . If the peak epochs are shifted by (r) the middle 
of the 90 percent confidence level on t , the peak flux - peak 
epoch best-fit becomes F p cc (t p — r) -2 ' . The probability 
for the F p — t p correlation to occur by chance decreases by a 
factor 3 when an afterglow time-shift (relative to the GRB 
trigger) is allowed. 

We have also identified 14 X-ray afterglows (Fig 4) that 
display a light-curve peak. Although XRT has monitored 
several hundreds of X-ray afterglows so far, the X-ray peaks 
are very rarely seen because they are, most often, overshined 
by the prompt GRB emission. The fc-corrected X-ray peaks 
also display a peak flux - peak epoch correlation: F p cc t~ 16 , 
which is slower than for optical peaks, and whose index 
changes little if a shift of the afterglow beginning is allowed. 
There is a 50 percent chance that the distributions of optical 
and X-ray peak epochs (Fig 5) are the same. 

The correlations identified for optical and X-ray peaks 
represent the bright and late edge of the entire distribution 
of peaks in the F p — t p plane, as peaks that are dimmer 
or occur earlier are more likely to be missed. Perhaps, the 
only model-independent conclusion that can be drawn from 
the existence of a bright and late edge in the distribution of 
peaks is the existence of an upper limit to the energy that 
afterglows radiate. Two thirds of optical and X-ray peaks 
decay faster than t^ 1 after the peak, thus F p t p is a good 

t Fortunately, allowing for a non-zero to does not change signifi- 
cantly the slopes of the optical and X-ray peak flux - peak epoch 
relations 



measure of the radiative output for most afterglows. Then, 
the peak correlation having a slope larger than unity (1.6 in 
X-rays, 2.0 in the optical) implies that later peaks radiate 
less energy than earlier ones. 

To reach more detailed conclusions, we have investi- 
gated two likely models for light-curve peaks, both pertain- 
ing to the synchrotron emission from the forward-shock the 
energizes the burst ambient medium. The slope of the bright- 
edge of the Fp — tp distribution enables a test of each model 
if it is assumed that that anticorrelation is induced by the 
variation of one parameter among the set of afterglows with 
peaks. 

One model is that of an off jet-aperture observer lo- 
cation, the afterglow peak occurring when the jet has de- 
celerated enough that its emission is relativistically beamed 
toward the observer. In this model, jets seen at a larger angle 
peak later and dimmer. We have derived the F p — t p correla- 
tion induced by a variation of the observer offset angle, for 
both conical and laterally-spreading jets, and for a homo- 
geneous or a wind-like medium, and we have found that a 
conical jet (and either type of medium) may explain the ob- 
served X-ray peak correlation, given the measured average 
X-ray spectral slope f} x — 1.0. If the optical peaks correla- 
tion is, indeed, steeper than for the X-ray peaks, that would 
require a wind-like ambient medium. 

The other model for afterglow peaks considered here is 
the onset of the forward-shock deceleration, occurring when 
the reverse shock has crossed the ejecta shell, the afterglow 
peak being due to a change in the dynamics of the forward- 
shock, caused by a decrease in the rate at which energy is 
transferred from the ejecta to the forward-shock. We have 
derived the F p —t p correlation expected for a semi-relativistic 
reverse-shock (F p and t p depend on the ejecta initial Lorentz 
factor and the ambient medium density) and for a relativistic 
reverse-shock (F p and t p depend only on the geometrical 
thickness of the ejecta shell) and have found that this model 
can explain the F p — t p relation measured for X-ray peaks if 
the ambient medium is homogeneous and if X-ray is below 
the synchrotron cooling frequency (for either type of reverse- 
shock), but it underpredicts the slope of the optical peak flux 
- peak epoch relation. 

There are three important caveats for any conclusion 
drawn from the slope of the optical and X-ray afterglow 
peaks. 

The first caveat is that the sample of X-ray peaks is 
small, hence the slope y x = 1.6 of the F p — t p relation is 
still quite uncertain. The X-ray afterglow 060614 displays a 
plateau around the epoch of the optical peak and, thus, was 
not included in the sample of X-ray peaks. However, if an 
X-ray peak is assigned at the epoch of the optical peak, then 
the best-fit to the peak fluxes and epochs of the resulting 
set of 15 X-ray afterglows would have a slope j x = 2.1, 
compatible with that measured in the optical (p/ — 2.0). 

The second caveat is that requiring a model to account 
for the slopes of both optical and X-ray peaks is justified 
only if the afterglow emissions at these two frequencies arise 
from the same population of electrons. A simple test for the 
common origin of the optical and X-ray afterglow emission 
is the simultaneity of the light-curve peaks. Most of optical 
peaks occurred while the X-ray afterglow was dominated 



by the prompt emission, hence the achromaticity of optical 
peaks is, in general, impossible to establish. The second op- 
tical peak of 100901 appears also in the X-ray, the optical 
and X-ray peaks of 100418 seem achromatic, but those of 
110213 are chromatic. There are four optical peaks (060206, 
060614, 070802, 081029) for which the X-ray data suggest 
the existence of a simultaneous peak, although they do not 
show it clearly, and one (110205) achromatic optical peak. 
For two X-ray peaks (080310, 080319C), the optical data al- 
low for a simultaneous optical peak occurring at the end of a 
light-curve plateau. All in all, there is evidence that optical 
and X-ray light-curve peaks occur simultaneously more of- 
ten than not, i.e. in favour of a common origin of the optical 
and X-ray afterglow emissions. 

The third caveat is that we have tested only the ability 
of one model parameter to induce the observed F p — t p an- 
ticorrelation but, in each model, there are other parameters 
that could alter the slope of that correlation. For the off- 
aperture observer model, where the fundamental parameter 
that drives the peak correlation is the observer's offset an- 
gle, the jet energy per solid angle and the ambient medium 
density may be seen as secondary parameters whose varia- 
tion among afterglows only introduces some scatter around 
the peak correlation induced by the main parameter. How- 
ever, for the deceleration-onset model, there are two basic 
parameters that can induce an F p — t p anticorrelation (the 
ambient density and the ejecta initial Lorentz factor, or the 
ejecta shell thickness), the resulting slope of the F p —t p rela- 
tion depending on the (unknown) width of the distribution 
of those two parameters. 

Moreover, given that both models for afterglow peaks 
considered here yield F p — t p relation slopes that are either 
smaller or larger than measured for X-ray and optical peaks, 
it is possible that a peak relation of intermediate slope re- 
sults from combining both peak mechanisms, some peaks 
being due to an off-jet observer location while others are 
caused by the onset of deceleration. That would make it 
even harder to use the F p — t p relations identified here for 
model testing. 
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